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Abstract 

In this talk, I explore the possibilities of probing the multiparton 
correlation functions at CEBAF at its current energy and the energies 
with its future upgrades. 

I. INTRODUCTION 

The perturbation theory of Quantum Chromodynamics (QCD) has been 
very successful in interpreting the data from high energy experiments. For 
observables involving hadrons, QCD factorization theorem enables us to 
systematically separates the short-distance partonic dynamics from the 
long-distance non-perturbative hadronic informations. In particular, when 
the energy exchange in a collision is high enough, the QCD factorization 
theorem factorizes all non-perturbative information into a set of universal 
parton distributions. Such universal parton distributions have been suc¬ 
cessfully extracted from the existing data. 

However, the nonperturbative dynamics of the strong interactions is 
much richer than what we have learned through the parton distributions. 
The parton distributions are interpreted as the probability densities to find 
the partons within a hadron. Knowing only such probability densities is 
not enough to understand the full partonic dynamics within a hadron. But, 
very little has been learned about the multiparton correlations. This is be¬ 
cause the physics associated with the multiparton correlations is directly 
related to multiple scattering, and most physical observables measured in 
existing high energy experiments are dominated by the contributions from 
the single hard scattering. And, in general, the multiple scattering is power 
suppressed in comparison with the single scattering. Therefore, it is often 
the case that even an excellent data cannot provide a good enough informa¬ 
tion for extracting the multiparton correlations, after removing the large 
contributions from the single scattering. 
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In this talk, I use examples to show that (1) by taking the advantage of 
kinematics, we can enhance the contribution from the multiparton scatter¬ 
ing; and (2) there are new types of physical observables which eliminate the 
contribution from the single scattering. The rest of this talk is organized as 
follows. A general discussion on relationships between the multiple scatter¬ 
ing, the multiparton correlation functions and the high twist operators is 
given in next section. In Sec. Ill, with some simple assumptions, we show 
analytically that when xb is large, the complete twist-4 contributions to 
the DIS structure functions can be simplified and is proportional to the 
derivative of the twist-2 parton distributions. In Sec. IV, we demonstrate 
how to extract multiparton correlation functions inside a large nucleus. Fi¬ 
nally, in Sec. V, we provide our summary and discussion on why CEBAF 
is a good place to measure the multiparton correlation functions. 

II. MULTIPARTON CORRELATIONS AND HIGH TWIST 

Multiparton correlation functions are defined as the Fourier transform 
of the matrix elements of the multiparton fields, 

r n ^ t 4 

cx / II WY elkm ' Vm 

x (^ r l0(O)^(2/i)---^(yn-i)|A r ), (1) 

where | N) is the hadron state, the ki with i = l,..,n — 1 are independent 
parton momenta, and the efts are the parton field operators. Because of 
the momentum dependence of the correlation functions, the multiparton 
field operators in the position space, ... <f>(y n -i) are not local. 

In a gauge theory, like QCD, a path integral of the gauge field is often 
necessary to be sandwiched between the </>’s in Eq. ([!]), in order to make the 
correlation functions gauge invariant jij . The simplest correlation functions 
are the matrix elements of operators with two parton fields. For example, 
the quark distribution of a proton of momentum p , moving in the “+” 
direction, 

q{x) = / (^) Mv~)\p) , (2) 

which is a correlation of two quark fields. 
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Beyond two-parton correlation functions, the choice of the multiparton 
field operators is not unique because of the spin structure of the quark 
fields, the equation of motion, and the gauge invariance. For example, the 
structure functions measured in deeply inelastic scattering (DIS) can have 
following factorized form, 

F(x b ,Q 2 ) = '^2c( m \x B ,x 1 ,x 2 ,..,Q‘ 2 /lJ?)®T( m \x 1 ,X2,..,iJ?) , 

m / 

(3) 

where x B is the Bjorken variable, and Q 2 is the large virtual momentum 
exchange in DIS. In Eq. (Q), the coefficient functions, C^’s, are calculable 
within the QCD perturbation theory, and the T^’s are nonperturbative 
and defined as the Fourier transform of the matrix elements of the multipar¬ 
ton field operators. The term with m = 0 is often called the leading twist 
or leading power contribution, while the terms with m 0 are called the 
high twist contributions or power corrections. It was pointed out in Ref. j^] 
that at a given order of power corrections, m, the T can be expressed 
in terms of the multiparton correlation functions with the same number of 
physical parton fields (i.e., fields with physical polarizations). For example, 
the T ^ for the m = 1 power corrections in DIS can be expressed in terms 
of four -parton correlation functions only (note, the power corrections due 
to the hadron mass are not included in discussion here). 

Although the /crar-parton correlation functions correspond to the lead¬ 
ing power corrections to the structure functions in DIS or other physical 
observables, we have not had much success in extracting the information on 
these correlation functions. This is because of following obvious difficulties: 

• The twist-4 (or higher twist) contributions are power suppressed in 
comparison with the corresponding leading twist contributions. 

• There are many more four-parton correlation functions than the two- 
parton correlation functions (or parton distributions) due to the rich 
spin decompositions and flavor combinations of different parton field 
operators. 

• Dependence on extra parton momenta requires much more informa¬ 
tion to extract multiparton correlation functions than the parton dis¬ 
tributions. 
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In order to extract the multiparton correlation functions, we have to identify 
the physical observables such that above difficulties can be minimized. For 
example, we can measure a physical observable in a particular part of the 
phase space so that the power corrections are relatively enhanced; we can 
identify an observable which gets no leading twist contribution; we can take 
an advantage of kinematics and identify observables with a limited leading 
subprocesses so that only a small number of multiparton correlation func¬ 
tions are relevant || [|] . In the rest of this talk, we will give two examples 
to show how to extract the four-parton correlation functions, and provide 
the physical interpretation of these four-parton correlation functions. 

III. TWIST-4 CONTRIBUTION TO THE STRUCTURE 
IN DEEPLY INELASTIC SCATTERING 

In this section, we demonstrate how to extract the leading power cor¬ 
rections (i.e., the twist-4 contributions) to the structure functions in DIS 
at the CEBAF energies. In particular, we provide the physical interpreta¬ 
tion of the leading power corrections. Note, the inclusion of target mass 
corrections are straightforward, and will not be discussed here. 

Although the twist-4 contributions to the structure functions in DIS 
are power suppressed by a factor of (l/Q 2 ), as shown in Eq. (|3|), the corre¬ 
sponding coefficient functions C^(xb, Xi, x 2 , ■ ■ ■) are larger than the lead¬ 
ing C(°)(xb)’s for large Xb ■ The structure functions at large xb are useful 
for extracting the information on the size of power corrections (^]. How¬ 
ever, after removing the leading twist contributions, the existing data are 
still not good enough for extracting detailed information on different four- 
parton correlation functions T^’s in Eq. (|Ij). It is a common practice to 
use the following parameterization for the leading power corrections ||, 

F 2 {x B ,Q 2 ) ~ ^1 + -qT^ F 2 (x b ,Q 2 )lt , (4) 

where F 2 (xb,Q 2 ) lt includes the leading twist (or leading power) con¬ 
tribution and the target mass corrections to the full structure function 
F 2 {xb,Q 2 )- The h(x) in Eq. (j|) is a fitting function. Yang et al. || used 
following parameterization, 

h{x B ) = a (—^- - c\ , (5) 

\ 1 xb ) 
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to fit the existing DIS data and extract the values of the parameters a , b 
and c for different targets. On the other hand, a complete expression of the 
twist-4 contributions to the DIS structure functions at the leading order of 
a s were derived more than fifteen years ago Hi- It is important to find 
the linkage between the phenomenological fits and the theoretical formulas, 
and reveal the physical meanings of the these fitting parameters. 

In terms of one photon exchange, the invariant mass square of the 
hadron-photon system is given by 

W 2 = [p + q) 2 ~ — (1 - x B ) ■ (6) 

x B 

For inclusive DIS, it is obvious that the power expansion of the structure 
functions in Eq. (g) should have terms proportional to (1 /W 2 ) m . Because 
W 2 shown in Eq. (ij) can be much smaller than Q 2 in large x B region, the 
power corrections are relatively more important in this region. Without 
worrying about the target mass corrections, the coefficient functions should 
have following asymptotic behavior, 

C^\x B , Xl ,x 2 ,...,Q 2 /^)=^ (—!—) (7) 

V1 — x B ) 


as x B — i > 1. From above simple argument based on the DIS kinematics, it 
seems natural to parameterize the fitting function h(x B ) in Eq. ©• 

In addition to the asymptotic behavior of the coefficient functions, it 
is also important to understand why the the twist-4 part of the structure 
function, ^(^b, Q 2 )tw 4 > which are proportional to the four-parton corre¬ 
lation functions, can be approximated to be proportional to the twist-2 part 
of the structure functions, as shown in Eq. (f|). The twist-4 contributions 
to the structure functions at the leading order of a s are given by !•§ 


F 2 (x b ,Q 2 ) TW4 = ^^2 e q 4T g { ^i(x s ) - x B J dxidx 2 

J(x,-x B )-S(x 1 -x {8) 


x 2 - X\ 


and 


F l {x b ,Q 2 ) tw 4 = 7« e 9 4 T qDi ( Xb ) 


Q 


(9) 
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plus the contributions from the four-quark correlation functions. In Eqs. (||) 
and (|]), the quark-gluon correlation functions are defined as 

T qD 1 ( x B) = 2 J ^—e iXBP+v (j>\ijj q (O)j a j + i 0 D%{O)D^(y-)i> q (y-)\p)-, 

( 10 ) 

and 

r(!) ( Xn ) — 1 [ dy P + dy I ix lP +y- i(x 2 -x 1 )p + y7 

Ud^b) - 4 J 2j: 2jt e e 

x(p|^g(o)7a7 + 7/3-D?(yr) £) T(yr)V’ g (y _ )b), (n) 


where Dt are transverse components of the covariant derivatives. 

In Ref. jij , we show that when xb is large, a simplified expression, simi¬ 
lar to the phenomenological parameterization, can be derived from Eq. (| 8 |). 
Key arguments and assumption of the derivation are summarized as follows: 

• It is argued that the correlation function 7’^ (xb) cannot be propor¬ 
tional to F 2 (xb,Q 2 )ut/(X ~ xb) as Xb —> 1- Therefore, the leading 
contribution to F 2 {xb, Q 2 )tw 4 is from the second term in Eq. (|j). 

• Assuming that D^(y^) = — g a y D^{y^)D^{y^) is a very slow func¬ 
tion of yi, we can reduce the quark-gluon correlation function T^ 2 
to the normal twist -2 quark distribution defined in Eq. (^J). 


• It is argued that the contributions from the four-quark correlation 
functions are smaller than that from the quark-gluon correlation func¬ 
tions given in Eq. (||). 

With the above arguments and assumption, we derive 


Fi(xb, Q 2 )tw4 




F 2 (xb,Q 2 ) TW 4 


xb 
2 (D* 


2Q 2 


xb 


- F L (xBiQ 2 )tW4: 


q{x B ), 


dXB 


( 12 ) 


where the twist-2 quark distribution q(xB) is defined in Eq. (|j). The {Dj) 
is an averaged value of the covariant derivative. Eq. ©) is a direct result of 
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above arguments and assumption. As shown below, Eq. ®D can be reduced 
to an expression very similar the phenomenological parameterization given 
in Eqs. (0) and (g). 

Eq. (|2|) in its present form is still different from the parameteriza¬ 
tion shown in Eq. (g) because of the derivative of the quark distribution. 
However, when x B is large, the quark distribution should have following 
behavior 

q{x B ) =>- (1 - x B ) aq asxB->l, (13) 

with the real parameter a q ~ 3.0 for the valence quark distributions. There¬ 
fore, we have 

x B — — q(x B )=a„ -—— q(x B ) + terms without -. (14) 

ax B l — x B l — x B 

Substituting Eq. (|b| ) into Eq. 

E 2 (a;_B, Q 2 )tw4 ~ 


It is now clear that our analytical result shown in Eq. ( |l2|) or equivalently 
Eq. (05h is very similar to the phenomenological parameterization shown in 
Eqs. (g) and (g). The difference is that our result was derived analytically 
from the complete leading order twist-4 contributions with some simple 
arguments and assumptions, and thus, every parameter have the physical 
interpretations. In addition, because of the dependence on the derivative 
of the quark distributions, instead of the leading twist structure function, 
our result provides the natural difference between the proton target and 
the neutron target, which is observed in the NMC data. 

In order to measure the power corrections to the structure functions, 
we like to keep x B large and Q 2 small. However, from Eq. (g), we cannot 
keep Q 2 too small, because of the hadronic resonances when W is less than 
2 GeV. When the invariant mass of the “hadron-photon” system is less 
than 2 GeV, the power expansion shown in Eq. (g) fails. Therefore, to test 
the perturbative power expansion of the structure functions and to measure 


(12), we obtain 


E 2 / \ a q{^r) x b 

e q x B q(x B )^^ - 


terms without 


2 Q 2 1 — x B 

1 


1 - 


(15) 
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the multiparton correlation functions, it is important to keep W > 2 GeV. 
With its luminosity and energy (in particular, with its future upgrades), 
CEBAF should be a good place to measure such power corrections, and to 
test the result shown in Eq. @- 

IV. MULTIPARTON CORRELATION IN A NUCLEUS 

Multiparton correlation functions inside a large nucleus are extremely 
important and useful for understanding nuclear dependence in relativistic 
heavy ion collisions. Inside a large nucleus, multiple scattering, which is 
directly associated with the multiparton correlation functions, can take 
place within one nucleon or between different nucleons. Since the leading 
single scattering at high energy is always localized within one nucleon, it 
cannot generate any large dependence on nuclear size. Similarly, multiple 
scattering within one nucleon cannot provide much dependence on nuclear 
size either. Therefore, large dependence on nuclear size is an unique signal 
of multiple scattering between nucleons. Measurement of such anomalous 
dependence on the nuclear size for any physical observable will provide 
direct information on the multiparton correlation functions in a nucleus. 
However, in order to test QCD, and to extract the useful information on the 
multiparton correlation functions, it is important to identify the physical 
observables which depend only on a small number of correlation functions. 
Otherwise, the data will not be able to separate contributions from different 
multiparton correlation functions. 

It was pointed out recently in Ref. J|] that at the leading order of a s , it 
needs only one type of multiparton correlation function for both transverse 
momentum broadening of Drell-Yan pair and the jet broadening in DIS. It 
is a quark-gluon correlation function, and is defined as 


t£ f {x b ) = 


X 


j d £e‘-‘^-&£ e(a -- yr)ll{ - y ; ) 

<C4in(»P«(o)^(r)f + “(!/r)ie4}, 


(16) 


with A the atomic weight. In terms of this correlation function, the jet 
broadening in DIS can be expressed as § 


A (P 2 T ) = 


47r 2 q s E g e q T£?( x b) 

3 E q ^q A (x B ) 


(17) 
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where pt is the transverse jet momentum in the photon-hadron frame in 
DIS. Measuring the jet momentum broadening, A(4), provides a direct 
measurement of the quark-gluon correlation functions inside a nucleus. 

However, at the CEBAF energies, it is not possible to measure the 
jet cross section and corresponding jet momentum broadening in nuclear 
targets. But, it should be possible to measure the transverse momentum 
broadening of the leading pions (e.g., final-state pions with a relatively 
large momentum). Define the averaged transverse momentum of the pions 
of momentum i as 


<4> eA 


f<u 2 T i 2 T 


d(JeA—* 7r 

dxsdQ 2 


ddeA 

dxsdQ' 2 


(18) 


and define the nuclear broadening as 


A(4) = (4) eA - (4) eJV . 


We obtain |J that 


A<4> = 


47 x 2 a a Eg eg i 1 Zmiri dzz 2 D q ^(z)Tf F (x B ) 
3 Eg e 2 q q A {x B ) 


(19) 


( 20 ) 


where D q ^ v {z) are the quark-to-pion fragmentation functions, and z m i n 
depends on minimum pion momentum and is frame dependent. Since the 
quark-to-pion fragmentation functions are known, Eq. ( |2C| ) provides the 
direct information on the quark-gluon correlation functions inside a nucleus. 
By measuring the transverse momentum broadening for both tt 21 and 7r°, 
and keep a reasonable large value of z min , we can extract the quark flavor 
dependence of the correlation functions. 


V. CONCLUSIONS AND OUTLOOK 

We present an analytical derivation of the leading twist-4 contributions 
to the DIS structure functions when xb is large, and argue that CEBAF is a 
good place to test such contributions. Our derivation provides a direct link¬ 
age between the complete twist-4 contributions and the phenomenological 
parameterization used for fitting the existing data; and physical meanings 
of the parameters in the phenomenological parameterization. Besides the 
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similarities between the analytical result and the parameterization, there 
are also some differences. It is the difference that makes the measurements 
at CEBAF even more interesting. 

We also demonstrated that measuring the transverse momentum broad¬ 
ening of the leading pion production at CEBAF can provide direct infor¬ 
mation on multiparton correlation functions in a large nucleus. Such infor¬ 
mation are extremely important for the multiple scattering and anomalous 
nuclear dependence in relativistic heavy ion collisions. 
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